ON TRIPLY EVEN BINARY CODES 



KOICHI BETSUMIYA AND AKIHIRO MUNEMASA 

Abstract. A triply even code is a binary linear code in which 
the weight of every codeword is divisible by 8. We show how two 
doubly even codes of lengths mi and TO2 can be combined to make 
a triply even code of length mi + m2 , and then prove that every 
maximal triply even code of length 48 can be obtained by combin- 
ing two doubly even codes of length 24 in a certain way. Using 
this result, we show that there are exactly 10 maximal triply even 
codes of length 48 up to equivalence. 



1. Introduction 

For the past few decades, extensive research of doubly even binary 
hnear codes has been done. These codes turned out to be connected 
with objects in various areas, for example, sphere packing problem, 
combinatorial designs, finite groups, integral lattices, modular forms 
and so on |H [18]. In this paper, we are concerned with a subclass of 
the class of doubly even codes, called triply even binary codes. A triply 
even code is a binary linear code in which every codeword has weight 
divisible by 8, in other words, a binary divisible code of level 3 in the 
sense of [13] . Dong, Griess and Hohn [5J pointed out that a certain 
triply even binary code of length 48 arose naturally from a Virasoro 
frame of the moonshine vertex operator algebra V'^. Subsequently, 
Miyamoto |T6] found a construction method of l^'' from that code. 
Lam and Yamauchi [12] formulated this construction for the class of 
framed vertex operator algebras. To be precise, a holomorphic framed 
vertex operator algebra of central charge n is constructed from a triply 
even code of length 2n whose dual is even. Unlike doubly even codes, 
the classification of all triply even codes of modest lengths has not been 
established yet. 

The purpose of this paper is to develop a basic theory of maximal 
triply even codes, and to give a classification of maximal triply even 
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codes of length 48. Since any triply even code of length up to 48 can be 
regarded as a subcode of some maximal triply even codes of length 48, 
one can derive easily the classification of all triply even codes of lengths 
up to 48. It turns out that every maximal triply even code of length n 
with n = (mod 8) and n < 40 is obtained as the generalized doubling 
D(C) of a maximal doubly even code C (see Definition [7]), and = 48 
is the smallest length with n = (mod 8) for which there exists a 
maximal triply even code not equivalent to 1){C) for any doubly even 
self-dual code C. The unique maximal triply even code C{Tiq) of length 
48 not equivalent to generalized doublings is obtained by augmenting 
the code C(Tio) of length 45 generated by the adjacency matrix of the 
triangular graph Tiq. 

A principal consequence of our classification is that every structure 
code of the moonshine vertex operator algebra lies in the generalized 
doubhng T){C) of a doubly even self-dual code C of length 24. This fol- 
lows from O Proof of Lemma 5.5], together with the fact that C(Tiq)-^ 
has minimum weight 2. Moreover, our result constituted a foothold for 
further developments toward the classification of holomorphic framed 
vertex operator algebras of central charge 24. The headway of the de- 
velopments is as follows. By Lam and Yamauchi [12], every triply even 
code of length a multiple of 16 containing the all-ones vector is the 
structure code of some holomorphic framed vertex operator algebra. 
So it is natural to ask whether there are any new holomorphic framed 
vertex operator algebras of central charge 24 having one of the codes 
we found as the structure code. Lam [10] recently constructed ten new 
holomorphic vertex operator algebras of central charge 24 using sub- 
codes of CiTio). Subsequently, Lam and Shimakura [TT] constructed 
seven new holomorphic vertex operator algebras of central charge 24 us- 
ing subcodes of D^diQ^e^) and D{ef^). They also obtained a complete 
classification of the Lie algebra structures of holomorphic framed vertex 
operator algebras of central charge 24. These holomorphic framed ver- 
tex operator algebras correspond to some of the conformal field theories 
predicted to exist by Schellekens [19j . 

This paper is organized as follows. In Section [21 properties and some 
construction methods of triply even codes are given. In Section 121 we 
prove that some maximal triply even codes can be constructed from 
doubly even self-dual codes by the doubling process. In Section HI an 
infinite series of maximal triply even codes is constructed by triangular 
graphs and some properties of the codes in this class are given. In 
Section [5l a method for constructing a triply even code from a pair 
of doubly even codes is given. The main result in Section [5] states 
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that every maximal triply even code is obtained from a pair of doubly 
even codes containing their radicals. In Section [HI an efficient method is 
described for determining whether a given doubly even code contains its 
radical. In Section [TJ we show that the method described in Section [5] 
gives all maximal triply even codes of length 48 and, as a result, a 
classification of maximal triply even codes of length 48 is given. In 
Section IHl a classification of maximal triply even codes of lengths 8, 16, 
24, 32 and 40 is given. The result is available electronically from [Ij, 
as well as a complete program in Magma |2] needed to produce it. 

2. Basic constructions for triply even codes 

Throughout the paper, a code will mean a binary linear code, or 
equivalently, a linear subspace of the vector space over the field F2 
of two elements. The support of a vector u = {ui, . . . , m„) G is the 
set supp(u) = {i \ Ui = 1}, and the weight of u is wt(M) = | supp(m)|. A 
triply even code is a code in which every codeword has weight divisible 
by 8. A doubly even code is a code in which every codeword has weight 
divisible by 4. 

In this section, we give basic properties of triply even codes, and 
construction methods of triply even codes from doubly even codes. An 
[n, k] code is a code C C Fg with dimC = k, and n is called the length 
of C. For codes C and D of length n, C is equivalent to D ii C = D" 
for some coordinate permutation a G S'„. The automorphism group 
Aut(C) of C is defined as {cr G S'„ | C = C^}. The linear span of a 
subset S* C F2 over F2 is denoted by {S). For u,v E Fg, we define 
M * w to be the vector in Fg with supp('U *v)= supp(m) fl supp(f ). For 
C,D C ¥2, we define C * D := {u * v \ u E C,v E D). For vectors 
u E F^ and v E Fj, we denote by (n | t>) G F™"*"" the vector obtained 
by concatenating u and v. For subsets C C Fj^ D C F2^, we define 
the direct sum of C and D as 

C®D = {{u\v)E F^^+"2 \ uEC,vED}. 

If C and C (resp. D and D') are codes of length ni (resp. 722) then 
{C ®D)* {C © D') = {C* C) ®{D* D'). A code C is said to be 
decomposable if it is a direct sum of two codes. We denote by G Fg 
and On G F2, the all-ones vector, the zero vector, respectively. We will 
omit the subscript if there is no confusion. 

For vectors u = {ui, . . . , m„), f = (fi, . . . , f„) G Fg, we denote hj u-v 
the standard inner product XliLi '^^^ dual code of a code C is 
defined as {m G Fg | m ■ f = for any v E C} and is denoted by C^. A 
code C is self-dual (resp. self- orthogonal) if C = C"*- (resp. C C C^). 
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There exists a doubly even self-dual code of length n, if and only if n 
is divisible by 8. If C and D are codes, then (C © D)^ = ® D^. 

The following lemma is a special case of [21t Theorem 5.3] (see also 
[m Proposition 2.1]). 

Lemma 1. Let C = (S) be a code generated by a set S . Then C is 
a triply even code if and only if the following conditions hold for any 
u,v,w & S: 

(1) wt{u) = (mod 8), 

(2) wt{u*v) = (mod 4), 

(3) wt{u*v*w) = (mod 2). 

Definition 2. Let C be a doubly even code of length n. We define 
functions 

Q:C — ^ F2, ^^^Y^ mod 2, 

B:CxC^^¥,, {v,u) ^ "^^^Y""^ mod 2, 

T : X X — y F2, (m, v, w) ^ wt{u *v*w) mod 2. 
Clearly, the following equalities hold: 

(4) Q(x + y) = Qix) + Qiy) + Bix,y) ix,yEC), 
(5) 

B{x, y + z) = B{x, y) + B{x, z) + T{x, y, z) {xeC, y,ze C^), 

(6) 

B{x + y,z) = B{x, z) + B{y, z) + r(x, y, z) {x,y G C, z e C"*"), 

(7) T{x,y,z) = {x,yeC,ze{C*C)^). 

The doubly even radical radC, and the triply even radical RadC are 
defined as 

radC = {yeC^\ B{x,y) = (Vx G C)}, 
RadC = e radC I Q{x) = 0}. 

Clearly 

(8) rad(C © D) = rad C © rad D, 

(9) Rad(C © L>) D Rad C © Rad D 
hold. 

In general, the radicals radC, RadC are not linear and not neces- 
sarily contained in C, even if C is triply even. An example is C = 
(Ig) © (Is)- However, the following holds. 
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Lemma 3. Let C be a doubly even code. Then radC C (C * C)"*". 

Proof. We note that (C * C)-*- = {z E \ T(x, y, z) = for any x,y E 
C}. Suppose x,y E C and z G radC. Since x + y E C, we have 
T{x,y, z) = by ([6]). Thus z E {C * C)^, and the result follows. □ 

Lemma 4. Let C be a doubly even code, and suppose x,y E radC. 

(i) If y E C , then x + y E radC. 

(ii) If X + y E C , then x + y E rad C. 

Proof. Observe that, by Lemma [3], x E {C * C)-^ holds. For any z E C, 
we have B{z,x + y) = T{x,y,z) by ([5]). If y E C, then T{x,y,z) = 



0. Thus (i) holds. If x + y E C, then T{x,y,z) = T{x,x + y,z) + 



T(x, X, z) = 0. Thus (ii) holds. □ 



Lemma 5. Let C be a doubly even code, and suppose x E radC and 
z E RadC. Then 

(10) a; + CnradC= (a; + Cn(C*C)^)nradC, 

(11) z + CnRadC ={z + Cn{C*C)^)r]RadC. 

Proof The containment x + C fl radC C (x + C fl (C * C)-^) fl radC 
follows from Lemma [3] and Lemma [^i)[ As for the reverse containment, 
suppose y E C r\{C * C)^ and x + y E rad C. Since x E rad C, we have 
y E rad C by Lemma l ^ii) Thus x + y E x + Cfl rad C and ffTOj) holds. 
From (flOl), 



{z + Cn{C* C)^) n RadC = (^ + C n radC) n RadC. 

Suppose y E C n radC. Since wt(z) = (mod 8), z + y E RadC if 
and only if wt(t/) = (mod 8). Therefore 

(^ + C n (C * C)^) n RadC = ;z + C n RadC. 



Thus (ii) holds. □ 



Lemma 6. Let C be a doubly even code and D = {C * C)^ fl C. Then 
the restriction B\cxD of B to C x D is a bilinear pairing and QId 
is a quadratic form with associated bilinear form i?|_DxD- Moreover, 
C n radC and C H RadC are linear subcodes of C . In particular, if 
radC C C (resp. RadC EC), then radC (resp. RadCj is linear. 

Proof. First, note that since C C C * C, we have D G {C * C)^ C C^. 
For any x,y E C and 2; G -D, we have T{x,y,z) = by ([7]), hence 
B{x+y, z) = B{x, z)+B{y, z) by (|6]). Also, for any x E C and y,z E D, 
we have T(x, y,z) = by (^^, hence B{x, y + z) = B{x, y) + B{x, z) by 
(|5]). Therefore, S is a bilinear pairing on C x D, and Q Id is a quadratic 
form with associated bilinear form -Bji^xD by (jl]). 
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By Lemma [3], C fl radC = {y E D \ B{x,y) = for any x G C}. 
Since B\cxd is linear in the second variable, C fl radC is a linear 
subcode of C. 

Also, by dH), Q is linear on C n radC. Then, C fl RadC = {x E 
C n radC | Q{x) = 0} is a linear subcode of C. □ 

Definition 7. Let C be a code of length n and set R = C (1 RadC. 
We define the extended doubling D{C) and the generalized doubling 
I)(C) as 

(12) 2)(C) = ((1„|0„), (0„|1„), {ix\x) \xeC}), 

(13) ViC) = {R®On,{{x\x)\x eC}). 

We note that if C is a doubly even code, then D(C) is a triply even 
code and 

(14) dimi)(C) = dimC + dim(CnRadC). 

Note also that if C is a doubly even [n,d] code and n = (mod 8), 
then I'(C) is a triply even code of length 2n, dimension d+1 or d + 2, 
depending on 1 G C or not. In particular, if C is a doubly even self-dual 
code of length n, then D(C) is a triply even [2n, n + 1] code. This is a 
particularly important construction in connection with framed vertex 
operator algebras and lattices (see |^). In the next section, we give a 
sufficient condition for C under which 2)(C) is a maximal triply even 
code. 

3. Maximality of triply even codes 

In this section, we discuss maximal triply even codes, that is, triply 
even codes not contained in any larger triply even code. 

Lemma 8. If C is a triply even code, then C C RadC. Moreover, 
equality holds if and only if C is a maximal triply even code. 

Proof. The first part is immediate from Lemma [H For a vector x. 
Lemma [T] implies that (C, x) is a triply even code if and only if x G 
(C * C)-*" n RadC = RadC by Lemma [31 Thus the result follows. □ 

Lemma 9. Let C = ©^^^ Cj be a maximal doubly even code where Ci 
is an indecomposable component of length Ui for i = 1, . . . ,k. Then 

k 

(15) radC = 0(s,), 

1=1 

k 

(16) RadC = 0(t,). 
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where 




In particular, if C is a doubly even self-dual code, then 



k 



(17) radC = RadC = 0(l„J, 



i=l 

k k 



(18) i)(c)^0I)(a) = 02)(a). 



i=l i=l 



Proof. By (jS]), it suffices to prove (fT^ when C is indecomposable. Sup- 
pose V G rad C and x G C. Then 



ByLemma|3l v G {C*C)^. Then for any ?/ G C, = v-{x*y) = {y*x)-y. 
Hence 



By f lTOj) and f l20|) . (C, * a;) is a doubly even code. By maximality, 
f * a; G C. Also since x G C was arbitrary, C is the direct sum of codes 
supported by supp(f) and its complement. Since C is indecomposable, 
we obtain v G (1). Hence radC C (1). Therefore (fT5i) holds. 

We claim that there is at most one i such that = 4 (mod 8). If 
there are distinct j such that i = j = A (mod 8), then Ci © Cj is 
not a maximal doubly even code. This contradicts maximality of C. 
Therefore (fT6|) follows from (fT5|) . 

If C is a doubly even self-dual code, then each Ci is a doubly even 
self-dual code, hence rii is divisible by 8. Now, (IT7|) follows from (ITBl) 
and ([TB]). 

By (IT7I) . we have 



(19) 



wt(f * x) =0 (mod 4) 



(20) 



f * X G C . 



CnRadC = 0(l„J 



2=1 



= 0a nRadCi 



j=i 
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and hence 2)(0^=iC'j) = 0i=i^(C'i)- Since Cj is indecomposable, 
(HZ]) implies i)(Ci) = This proves □ 

Proposition 10. For any doubly even self-dual code C , {T){C)*D{C))'^ = 
D{C). In particular 2)(C) is a maximal triply even code. 

Proof. Suppose that C is an indecomposable doubly even self-dual code 
of length 2n. Then ([18]) implies T>{C) * I)(C) = D{C) * 'D{C) = 
C©C + D(C*C), hence dim(I)(C) *I)(C)) = 3n-l = 4n-dim:D(C). 
This implies that * I)(C))^ = By ([18]), the identity 

holds also for any decomposable double even self-dual code C. Now 
rad2)(C) C 2)(C) by Lemma [21 and hence is a maximal triply 

even code by Lemma [HI □ 

Example 11. It is known that the [8, 4, 4] Hamming code eg = 2)((l4)"'") 
is the unique doubly even self-dual code of length 8, up to equivalence. 
Also, d^Q = ©((Is)"*") and es © es are the only doubly even self-dual 
codes of length 16, up to equivalence. By Proposition[TOl ©(eg), D((i^g), 
D{es®es) are maximal triply even code of dimension 5, 9 and 10 respec- 
tively. In particular ©(eg) = ©(eg) is the Reed-Muller code RM(1,4) 
and T){es®es) = RM(1,4)®2. 

Example 12. It is known [I7j that there are precisely 9 doubly even 
self-dual codes of length 24. Two of these 9 codes are decomposable, 
and they are df^ © eg and ef^. The remaining 7 codes are indecompos- 
able and they are denoted by g24, d24, df^ , {dioc'^)'^ , dl'^ , dp~ , d^'^ . By 
Proposition [T0[ D{C) is a maximal triply even code for any of the 9 
doubly even self-dual codes C. We note from ([T7)) that Rad C G C and 
dim Rad C is the number of indecomposable components. Thus, for in- 
decomposable doubly even self-dual codes C of length 24, dim2)(C) = 
13 holds. Also, i)(d^6 © eg) = ©(©((Ig)^)) © RM(1,4) has dimension 
14, while D(ef^) = RM(1,4)®3 has dimension 15. 

Remark 13. As shown in Example [T2l the dimension of maximal triply 
even codes varies even if the length is fixed. The largest possible di- 
mension of triply even codes, however, has been determined in [22j . 
and the codes achieving the largest dimension have been determined 
in [H]. 

4. Triply even codes constructed from triangular graphs 

Let n be a positive integer with n > 4, and let f2 be a set of n 
elements. We denote by (2) the set of two-element subsets of The 
triangular graph T„ has the set of vertices (2), and two vertices 
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are adjacent whenever |q; n /3| = 1. It is known [8] that the graph T„ is 
a strongly regular graph with parameters 

(Tl\Tl — 1 ) \ 
^^-^,2(n-2),n-2,4j . 

Let An denote the adjacency matrix of T„. Then every row of An has 
weight 2{n — 2), and for any two distinct rows of the size of the 
intersection of their supports is either n — 2 or 4. Let C(T„) be the 
binary code with generator matrix An- 

It is clear that the code C{Tn) is triply even only if n = 2 (mod 4). 
The converse also holds by the following lemma. 

Lemma 14 (Haemers, Peeters and van Rijckevorsel [HI Subsection 4.1]). 
Ifn = 2 (mod 4), the weight enumerator of C(Tn) is 

L(«-i)/4j . s 

In particular, C{Tn) is a triply even of dimension n — 2. 

Let ai = {i, n} G (2) , and we denote by rj the row of An indexed 
by ai i.e., {k,l} G supp(rj) if and only if \ai fl {A;,Z}| = 1. Then the 
following lemma holds. 

Lemma 15 (Key, Moori and Rodrigues P, Lemma 3.5]). Ifn is even, 
then {vi \ i = 1,2, . . . ,n — 2} is a basis of C{Tn) ■ 

We note that the dimension of C(T'„) has already been determined by 
Tonchev ^20^ Lemma 3.6.6] and Brouwer and Van Eijl [3j. An explicit 
basis of C(Tn) is needed in the sequel to establish maximality of C(T„). 
The weight enumerator given in Lemma [T3] can also be derived from 
the basis. 

Lemma 16. If n is even, then {vi * rj \ 1 < i < j < n — 2} is a basis 
of C (Tn) * C (Tn) ■ In particular, 

dim{C{Tn) * C{Tn)) = 

Proof. Observe that, foT 1 < i < j < n, we have 

(21) supp(rj * rj) = {{i,j}} U {a^ \ I < k < n, i,j}. 

Suppose 

n-2 

Qri + Y Cijrj * Tj = 0, 

i=l l<i<j<™-2 
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where Ci,Cij G F2. Then q = for i = 1, ... ,n — 2, because \ai fl 
{j, n — 1}| = 1 if and only if i = j. Thus 

Yl ^^'^^^ * = 0- 

i<j<i<"-2 

For i,j, k, I G {1, . . . , n — 2} with i ^ j, k ^ I, fl2Tl) implies {fc, /} G 
supp(rj * Tj) if and only if {k, 1} = {i,j}- This implies Cjj = 0. □ 

Lemma 17. Ifn = 2 (mod 4), then {C{Tn) * C{Tn))^ = C(Tn) + (1). 
In particular, C{Tn) is a maximal triply even code. 

Proof. By (C(T„) * C(T„))^ D C(T„) + (1) and comparing the di- 
mensions using Lemmas [15] and [161 we obtain (C(T„) * C(T„))-'- = 
C(T„) + (1). Since wt(l) = = 1 (mod 2), Lemma [3] implies 

Rad C(Tn) C C(Tn). Thus C(Tn) is a maximal triply even code by 
Lemma [H] □ 

We define C'(T!„) to be the code of length / = 8[| "'^"~^'' l constructed 
from C(Tn) together with the all-ones vector of length /, i.e., C(Tn) = 
{h) + C(T„,) © 0;, where /' = / - 

Theorem 18. Ifn = 2 (mod 4), then C{Tn) is a maximal triply even 
code. 

Proof. Let / = 8[|^^^]. Then 

(^(T^) * C{T^))^ = {{h) + (C(T„) * C(T„)) © 0)^ 



(i,)^n((C(T„)*c(r„))^©F^; 



= {h)^ n ((C(T„) + (1)) ©F^') (by Lemma [IZD 

= C(T„)©(l,)^ + (li) 

= ^(T^) + © (Ip)^. 

Since V < 8, Lemma [3] implies Rad C(T„) C C(T„). The result follows 
from Lemma [HI □ 

5. Triply even codes constructed from pairs of doubly 
even codes with isometries 

In Section [21 we gave construction methods for a triply even code 
from a doubly even code. In this section, we give a generalization of 
these construction methods for a pair of doubly even codes. 

For a set of coordinates {ii, ^2, • • • , it] C {1,2,..., n}, let vr : Fg — )■ 
F2, tt' : F2 — >■ F2~* be the projection to the set of coordinates {ii, ^2, • • • , it], 

{ Jl ) • • • ) in—t 

}, respectively, where {ji, . . .Jn-t] = {1, • • • , . . . , it}. 
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For a code C of length n, the punctured code and the shortened code 
of C on a set of coordinates {^i, ■ ■ ,it} are the codes vr'(C), {7r'(c) | 
c G C, 7r(c) = 0}, respectively. 

Let Ci and C2 be doubly even codes and Ri be a subcode of Ci fl 
RadCj for i = 1,2. A bijective linear map 

(22) / : Ci/i?i C2/R2 

is called an isometry if wt(xi) = wt(x2) (mod 8) for any Xi + Ri G 
Ci/Ri and X2 + R2 G /(a;i + -R1). We note that if x + Ri = y + Ri with 
x,y G Ci, then wt(x) = wt(?/) (mod 8). The set of isometries (p2l) is 
denoted by $(Ci//2i, Ca/i^z)- 

For an isometry / G $(Ci/-Ri, CaZ-Ra), "we define a code 

(23) D(Ci,C2,i?i,i?2,/) = {(a^ilxs) I xi G Ci, 0:2 G /(xi + i?i)}. 
Since / is a bijective linear map, 

(24) D(Ci,C2,i?l,i?2,/) = {{X1\X2) I X2 G C2, Xi G /-^(X2 + i?2)}. 

Proposition 19. Let Ci be a doubly even code of length rrii for i = 1,2 
and Ri be a subcode of Ci fl RadCj. If f & $(Ci/i?i, C2/-R2); then 
the code D{Ci, C2, -Ri, -R2, /) is a triply even code of length nii + m2 of 
dimension dim Ci + dim R2 = dim Ri + dim C2 ■ 

Proof. Fix Ci, C2, Ri, R2 and /. We abbreviate D{Ci,C2, Ri, R2, f) 
as D. Since / is a linear map, D is linear. Since Ci and C2 are doubly 
even codes and / is an isometry, all the weights of elements of D are 
multiple of 8, that is, D is a triply even code. Moreover 

\D{C,,C2,Ri,R2,f)\ = \Ci\ X \f{R,)\ = \C^\ x |i?2|. 

Therefore dimD = dimCi + dimi?2 = dimi?i + dimC2. □ 

Remark that the construction method in Proposition [19] contains the 
constructions D(C) in (fT2|) and D(C) in (fT3!) as special cases. Indeed, 
let C be a doubly even code of length n. Then we have 

(25) D(C) =D{C,C,Cn Rad C,Cn Rad C, id) 
and if, moreover, n = (mod 8), then 

2)(C) = D(C+(l),C+(l),(l),(l),id). 

Note that, given doubly even codes Ci, C2 and subcodes Ri C CiH 
RadCi, i?2 C C2 n RadC2, the set $(Ci/i?i, C2/-R2) may be empty, 
and in this case Proposition [19] produces no triply even codes. We shall 
give a necessary and sufficient condition for the set <I>(Ci/i?i, C2/-R2) 
to be non-empty in Proposition [21] below. First we need to introduce 
some terminology. 
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Definition 20. Let C be a doubly even code, and let i? be a subcode 
of CnRadC. Let 

X = {x + i? e C/i? I wt(x) = (mods)}. 

We call the elements of the set X singular points of C/R. We denote 
by Qi{C,R) the setwise stabilizer of X in GL{C/R). The triply even 
check code C{C, R) of (C, R) is defined as 

C{C, R) = {c= {c^ e¥2\x e X) \ c^x e R}. 

By the definition, Qi{C,R) acts on C{C,R) as automorphisms, but 
the action is not necessarily faithful. Indeed, X may not span C/R. 

Proposition 21. Let Ci he a doubly even code for i = 1,2 and Ri be 
a subcode o/Cj flRadCj. Suppose that dimCi/i?i = dimC2/-R2- Then 
C{Ci, Ri) = C{C2, R2) if and only if ^{C^/Ri, C2/R2) 7^ 0- 

Proof. If C{Ci,Ri) = C{C2,R2), then there exists a bijection / from 
the set Xi of singular points of Ci/Ri to the set X2 of those of C2/R2 
which induces an equivalence from C(Ci, to C(C2, -R2)- It follows 
from the definition of the triply even check code that the bijection / 
extends to a linear mapping (Xi) — {X2). Extending further to Ci/Ri 
in an arbitrary manner, we obtain an isometry from Ci/Ri to C2/R2- 
The proof of the converse is immediate. □ 

The next proposition shows that every triply even code can be con- 
structed by means of the construction described in Proposition [121 

Proposition 22. Let D be a triply even code of length n. Fix a code- 
word X E D of weight rrii with < rrii < n. Let Si = supp(l + x) 
and S2 = supp(x) and let Ci and Ri be the punctured code and the 
shortened code of D on Si, respectively, fori = 1,2. Then Ci is doubly 
even, Ri G CiH RadCj for i = 1,2, and 

D^D{Ci,C2,Ri,R2,f) 

for some f G $(Ci//?i, C2/i?2). 

Moreover, if D is maximal, then RadCj = Ri for i = 1,2. 

Proof. All the statement except on the last one follows easily from 
Lemma [H Let tti : ^ 712 : ^ F2""'' be the projection 

to the set of coordinates supp(x), supp(l + x), respectively. Define 
TT : ^ F^ by 7t{x) = (7ri(x)|7r2(x)) (x G F^). Then d = Tii{D) 
{i = 1, 2) and D is equivalent to 7r(D). It is clear that the mapping 

/ : Ci/Ri C2/R2 

ci + Ri {x e C2 I (ci|x) G n{D)} 
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is a well-defined isometry and tt{D) = D{Ci, C2, Ri, R2, /)• 

If D is a maximal triply even code, then so is Tr{D). This implies 
(ri|0), (0|r2) G Rad7r(L)) for ri G RadCi and G RadC2- By 
Lemma [HI Rad7r(D) C vr(D). Therefore Ri C RadCj for i = 1,2. 
Hence the result follows. □ 

Proposition [22] indicates that every triply even code of length n con- 
taining a codeword of weight mi can be constructed from a pair of 
doubly even codes of lengths mi and n — mi. We will classify maximal 
triply even codes of length 48 by setting mi = 24 in Section [71 

For fixed codes Ci, C2 and i?i C Ci n Rad Ci, C C2 n Rad C2 the 
resulting code 

-D(Ci, C2, -Ri, R2, f) 

depends on the choice of the isometry /. However, some of these codes 
are equivalent to each other. The first algorithm is to check this, that 
is, we will give a sufficient condition for two resulting codes to be 
equivalent. We need this algorithm to reduce the amount of calculation 
to be reasonable. 

First, we define some groups. For a code C and a subcode R C 
C n RadC, we denote by Qo{C,R) the subgroup of GL{C/R) induced 
by the action of Aut(C) fl Aut(i?) on C/R and denote by Gi{C, R) the 
subgroup ^{C/R,C/R) of GL(C/i?). By the definition, the group 
Go{C,R) is a subgroup of GiiC,R). If R = C D RadC, then we 
abbreviate Qo{C,R), Qi{C,R) as Qo{C), Qi{C), respectively. If / G 
$(Ci/i?i,C2/i?2), then 

(26) gi{C,,Ri) = f-'og,{C2,R2)of 
and 

(27) $(Ci/i?i, C2/R2) = fogi{C,, i?i) = 6^1 (C2, R2) o f. 

If we replace / by cr2 o / o en, where (jj G Qo{Ci, Ri), then the resulting 
codes are equivalent, that is, 

D(Ci, C2, R2, f) = D{Ci, C2, R2, ^2 o / o di). 

This means that, in order to enumerate 

{D(Ci, C2, Ru R2, h)\he $(Ci/i?i, C2/R2)} 

up to equivalence, we first fix / G <l'(Ci/i?i, C2/-R2), and it suffices to 
enumerate the codes D{Ci,C2, Ri, R2, f o g) where g runs through a 
set of representatives for the double cosets 

if-' o 6^0(^2, R2) o f)\giiCi, Ri)/go{Ci, i?i). 
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6. Doubly even codes containing their radicals 

In view of Propsition [221 it will be necessary to extract only those 
doubly even codes C which satisfy RadC C C, in order to enumerate 
maximal triply even codes. In this section, we will give a criteria to 
verify whether a doubly even code C contains its triply even radical 
i.e., RadC C C. 

Throughout this section, let C be a doubly even code containing 1, 
and we denote {C * C)-^ r\ C hy D . For x G C"*", one can define a 
mapping B.j. : C ^ ¥2 by B.j.{c) = B{c, x) (c G C). By B.j. is linear 
when X G (C * C)"*". Thus we obtain a map 

0: (C*C)^ ^ Hom(C,F2) 
X ^-J■ Bx- 

By Lemma [3l we can write 

(28) 0-^(0) =rad a 

We remark that the map is not linear in general. More precisely, if 
we define a bilinear map 5 as 

5 : X ^ Hom(C,F2) 

(x,?/) {v t-^ T{x,y,v)), 

then for x,y E {C * C)-*-, 

0(x + y) = (f){x) + (f){y) + 6{x, y) 
holds by ([n]). In particular, ([7]) implies 

(29) 0(x + y) = 0(x) + 0(1/) (xG (C*C)^, 2/GD), 

and is linear on D. 

The function Q from Definition [2] can also be defined on radC, so 
we denote it by the same Q as follows. 

Q:radC^F2, m mod 2. 

Then RadC = g-i(O), and 

(30) Q{x + y) = Q{x) + Q{y) (x G C n radC, y G radC). 
Lemma 23. For a coset M G {C*C)-^/D, the following are equivalent. 

(i) 0(M) n 0(D) ^ 0, 

(ii) 0(M) = 0(D), 

(iii) MnradC 7^ 0. 



Moreover, if C HradC 7^ C Pi RadC, then each o/(i) -(iii) is equivalent 
to 
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(iv) MnRadC ^ 



Proof. Equivalence of (i) -(iii) follows immediately from fl28|) and f l29|) . 



Suppose C n rad C C H Rad C. It suffices to show that (iii) implies 



;iv) 



Suppose X G M n radC. If Q{x) = 0, then clearly (iv) holds, so 
suppose Q{x) = 1. By assumption, there exists y & C flradC such 
that Q{y) = 1. Then x + y e M, (j){x + y) = 0(x) + (f){y) = by (129|), 



hence x + y E radC. Moreover, Q{x + y) = Q{x) + Q{y) = by fl5U]) . 



Thus X + y G RadC, and hence (iv) holds. □ 



Proposition 24. Let C be a doubly even code of length a multiple of 
eight, containing 1. Suppose CflradC 7^ CflRadC. Then RadC (t C 
if and only if there exists a cos et M G (C * C)^/D satisfying (j){M) n 
0(L>) andM i^D. 

Proof. Since RadC C radC C (C * C)^ by Lemma El RadC ^ C if 
and only if M fl RadC ^ for some coset M G (C * C)^/D different 
from D. The resut then follows from Lemma [231 □ 



In view of equivalence of (i) and (ii) in Lemma [23| one can check the 
condition 0(M) n0(D) 7^ by testing whether an arbitrarily chosen el- 
ement x G M satisfies (pix) G (piD). Thus, the above proposition gives 
a convenient criterion for Rad C C C in terms of coset representatives 
for (C * C)-^/D, provided C fl radC 7^ C fl RadC. In the case where 
C n rad C = C n Rad C, the situation is slightly more complicated. 

Lemma 25. Suppose C n radC = C n RadC, and M e {C * C)^/D. 
IfMnRadCy^ 0, then 

(31) MnradC = MnRadC. 

Proof. By assumption, there exists x G MflRadC. Then by Lemma O 
we have 

(32) MnradC = ,T + CnRadC. 

Since x G RadC, ([II]) implies x + C fl RadC C RadC, hence M fl 
rad C C Rad C by ([32]) . This proves M n rad C C M n Rad C, and the 
reverse containment is trivial. □ 

Proposition 26. Let C be a doubly even code of length a multiple of 
eight, containing 1. Suppose CflradC = CflRadC. Let {xi, . . . ,Xt} C 
{C*C)^ be a set of coset representatives for the cosets M E {C*C)-^ / D 
satisfying (j){M) fl 0(L') ^ and M ^ D. For each i G {1, . . . 
choose yi E D in such a way that 4>{xi) = (piyi). Then the following 
are equivalent. 
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(i) RadC t C, 

(ii) wt{xi + Hi) =Q (mod 8) for some i E {1, . . . , t}. 

Proof. First, we note that 4>{M)n(j){D) ^ implies that 0(M) = 0(D) 
by Lemma l23l Thus there exists Ui E D such that (j){xi) = (f){yi), no 
matter how we choose a representative Xi for the coset Xi + D. 



Suppose (i) holds. Take x G RadC \ C and set M = x + D. Then 
X E M n RadC, and hence fl3Tl) holds by Lemma [251 Also, as x ^ D 
and 0(x) = by fl28|) . M = Xj + Z) holds for some i E {1, . . . ,t}. Thus 
Xi + Hi G M, while 0(xi + yi) = 4>{xi) + (piyi) = by fl2^ . Therefore, 
a^i + 1/i £ ^ n radC C RadC by fl^ . This implies wt(xj + ^j) = 
(mod 8). 



Conversely, if (ii) holds, then Xi + G Rad C \ C, and hence 



1 



holds. □ 



7. Classification of maximal triply even codes of length 

48 

In this section, we aim to give a classification of maximal triply even 
codes of length 48. In Section [3] and Section HI we gave 10 distinct 
maximal triply even codes of length 48. Now we show that the list is 
complete for a classification up to equivalence applying Proposition [22] 
and [in] for n = 48 and mi = m2 = 24. To do this, we first need 
to establish the existence of a codeword of weight 24 in any maximal 
triply even code of length 48. 

Lemma 27. Let D be a maximal triply even code of length n. Let T 
be the graph with vertex set {1, . . . , n} and edge set 

{supp(x) I X G -D"*", wt(x) = 2}. 

Then the following hold: 

(i) every connected component of T is a complete graph with at 
most 8 vertices, 

(ii) if there is a connected component of T with more than 4 ver- 
tices, then any other connected component has at most 3 ver- 
tices. 

Proof. Since is a linear code, it is clear that every connected com- 
ponent of r is a complete graph. Suppose that there is a connected 
component i^' of F with \K\ > 8. Then there exists a vector x G 
with wt(x) = 8 and supp(x) C K. Since the restriction oi y to K is 
or 1 for any y E D, we have wt(x * y) = or 8. This implies that 
{D,x) is triply even. Taking i E supp(x) and j E K \ supp(x), the 
vector with support {i,j} belongs to and is not orthogonal to x. 
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Thus X ^ D. This contradicts the fact that D is maximal, and the 
proof of (i) is complete. 

To prove (ii), suppose that there are distinct connected components 
K,K' of r with \K\ > 4 and \K'\ > 4. Then there exists a vector 
X G with wt(x) = 8, | supp(x) fl i^l = | supp(x) H K'\ = 4. Since 
the restriction of y to or K' is or 1 for any y & D, v/e have 
wt(x * ?/) = 0,4 or 8. This implies that {D,x) is triply even. The rest 
of the proof is exactly the same as (i). □ 

Lemma 28. Let D be a maximal triply even code of length 48 contain- 
ing 1. Then D has at least one codeword of weight 24. 

Proof. By Lemma [271 the number of codewords of with weight 2 is 
of the form 




where the summation is taken over the set of connected components of 
the graph F defined in Lemma [271 Let Ai > A2 > ■ ■ ■ be the partition 
of 48 associated with the decomposition of the vertex set of F into 
connected components. Lemma [27| implies that one of the following 
holds: 

(i) 4 < Ai < 8 and Ai < 3 for all i > 2, 

(ii) Ai < 4 for all i > 1. 

It is not difficult to show that the maximum value of (2') is (2) + 
13 g) = 67 for the case (i), and 12(5 = ''2 for the case (ii). Therefore, 
we conclude that has at most 72 codewords of weight 2. 

Now suppose that D has no codeword of weight 24, so that its weight 
enumerator is 

+ aX^°yS + (2^-1 - (1 + a))(X3Vi6 + X'^Y^'') + aX'Y^' + 

where k = dimD. It follows from the Mac Williams identities that the 
number of codewords of weight 2 in D-^ is 

3 ■ 2^-'=a + 104 + 2^^-*^ 

which is certainly greater than 72. This is a contradiction. □ 

In order to construct all maximal triply even codes of length 48 
by means of Proposition [19] and [22] for n = 48 and nii = m2 = 24, 
it suffices to consider the codes of length 24 satisfying Rad Ci C Ci 
as candidates for Ci and C2. This is because, if a resulting code 
D{Ci,C2, Ri, R2, f) is maximal, then Ri = RadQ for i = 1,2 as we 
mentioned in Proposition [221 
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We are now ready to describe our enumeration using Magma system 

As the first step, we enumerate all doubly even codes of length 24 
containing its triply even radical. Since there is a database of doubly 
even codes we could make use of it and extract only those which 
contain the triply even radical. However, since every doubly even code 
is equivalent to a subcode of the nine doubly even self-dual codes of 
length 24 [T7], we can find all the desired doubly even codes by suc- 
cessively taking subcodes of codimension one starting from the doubly 
even self-dual codes. This approach has an advantage that once we 
encounter a doubly even code C with RadC ^ C, then RadC" ^ C 
for any subcode C of C, so that it is no longer necessary to consider 
subcodes of C by Lemma [HI Table [1] gives the numbers of doubly even 
codes of length 24 containing its triply even radical with each given 
dimension and dimension of its triply even radical. 

Table 1 . The numbers of doubly even code C of length 
24 with RadC c C 



dim C \ dim Rad C 


1 


2 


3 


4 


5 


6 


12 


7 


1 


1 











11 


33 


6 


3 











10 


130 


19 


10 


1 








9 


308 


40 


23 


5 





1 


8 


363 


37 


25 


10 


1 


1 


7 


180 


16 


10 


11 


2 


1 


6 


27 


2 





4 


2 


1 


5 














1 






As the second step, we enumerate all resulting codes 
L'(Ci,C2,RadCi,RadC2,/o^) 

obtained from the all combinations of doubly even codes Ci, C2 above 
and a representative g E g ioi each double coset 

g e if-' og,{C2,R2)o f)\g,ic^,R,)/goic^,R,), 

where / is a fixed element of $(Ci/ RadCi, C2/ RadC2) by the proce- 
dure given in Proposition [19] 

We denote the set of doubly even codes of length 24 by 

^ = {fi'245 '^24' ^12 y ('^1067)^, dp', dp', dl'^, d'^Q © eg, ef ^}, 

in accordance with the notation of |T7]. From the combinations with 
Ci = C2, we obtain 1482 triply even codes. However, many codes of 
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them of the form ( 125|) turn out not to be maximal. This is because, if 
there is a doubly even code C such that C C C" and RadC = RadC", 
then 2)(C) C X)(C"). Therefore, we find that only 216 codes among 
the 1482 codes are possibly maximal. Then we use Lemma [8] to check 
maximality, and we are able to confirm that only 30 codes among them 
are maximal. Each of these 30 codes turns out to be equivalent to D(C) 
for some C G A. 

From the combinations with Ci ^ C2, we obtain 225 triply even 
codes, and 5 codes among them are maximal. One code is equivalent 
to C(Tio). The other codes are equivalent to a member of | C G 

A}. Therefore we obtain the following theorem. 

Theorem 29. Every maximal triply even code of length 48 is equivalent 
to for some C e A or C{Tio). 

The computer calculation needed to establish Theorem [22] was done 
by magma [2] under the environment using Intel® Core^M 2 Duo CPU 
T7500 @ 2.20GHz, and it took 650.240 seconds. 

8. Classification of maximal triply even codes of lengths 

8, 16, 24, 32 AND 40 

In this section, we give a classification of maximal triply even codes 
of lengths 8, 16, 24, 32 and 40 by using a shortening process from 
the results of maximal triply even codes of length 48 in the previous 
sections. 

It is easy to see that every maximal triply even code of length n is 
a shortened code of a maximal triply even code of length n + 1. From 
the list of maximal triply even codes of length 48, we can derive the list 
of all maximal triply even codes of shorter lengths by the shortening 
process. The shortened code of D(C) on one coordinate has an odd 
length, so it cannot be of the form D{C') for any C . However, for 
lengths divisible by 8, the following holds. 

Theorem 30. For = 4, 8, 12, 16 and 20, every maximal triply even 
code of length 2n is of the form i>{C) for some maximal doubly even 
code C of length n. 

Table Ogives the numbers of the maximal triply even codes of lengths 
8, 16, 24, 32 and 40, up to equivalence. 

In Table [2], the first and fifth columns indicate the length of each 
doubly even code and each triply even code, respectively. The second 
and sixth columns indicate the dimension as well. The third column in- 
dicates the number of indecomposable components of the doubly even 
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code. The fourth and seventh columns indicate the number of codes 
satisfying the condition. The eighth column gives the other construc- 
tion method to obtain it. 

Note that if C is some maximal doubly even code and k is the num- 
ber of self-dual indecomposable components of C, then dimD(C) = 
dimC + k hy (fill) and (ITB]) . For example, there is a unique doubly 
even [20, 9] code C which is the direct sum of three indecomposable 
codes, two of which are self-dual. Then D{C) is a triply even [40, 11] 
code. Similarly, there is a unique doubly even [24, 12] code C which is 
the direct sum of three indecomposable self-dual codes. Then 2)(C) is 
a triply even [48, 15] code. 

Table 2. The numbers of maximal triply even codes of 
lengths multiple of 8 up to 48 



maximal doubly even codes 


maximal triply even codes 


len 


dim 


# compos 


#codes 


len 


dim 


#codes 


remark 


4 


1 


1 


1 


8 


1 


1 




8 


4 


1 


1 


16 


5 


1 


C{T,) 


12 


5 


1 


1 


24 


5 


1 








2 


1 




6 


1 




16 


8 


1 


1 


32 


9 


1 








2 


1 




10 


1 




20 


9 


1 


7 


40 


9 


7 








2 


2 




10 


2 








3 


1 




11 


1 




24 


12 


1 


7 


48 


13 


7 








2 


1 




14 


1 








3 


1 




15 


1 












48 


9 


1 
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/* 

Appendix A. A Magma program for classification 

Enumeration of doubly even codes of length 24. This appendix 
gives Magma scripts to verify Theorem [29] in Section [7l 

It is known that there are precisely 9 doubly even self-dual codes of 
length 24 up to equivalence [IT]. The object desd24genmats is the list 
of generator matrices in the hexadecimal expression. Also the object 
desd24 is the list of the codes 

A = {(/24, 44, dlt, (rfioe?)+, dl+, dt, dl+, dte © eg, ef }. 



1 desd24genniats: = [ 



2 


[ 0xC75001, 


0x49F002, 


0xD4B004, 


0x6E3008, 


0x9B3010, 


0xB66020, 


3 


0xECC040, 


0xlED080, 


OxSDAlOO, 


0x7B4200, 


0xBlD400, 


0xE3A800 ] , 


4 


[ 0x7FE801, 


0x802802, 


0x804804, 


0x808808, 


0x810810, 


0x820820, 


5 


0x840840, 


0x880880, 


0x900900, 


OxAOOAOO, 


OxCOOCOO, 


OxFFFOOO ] , 


6 


[ 0x7E0F81, 


0xFC0082, 


0xFC0104, 


0xFC0208, 


0xFC0410, 


0xFC0820, 


7 


0x820FC0, 


0x861000, 


0x8A2000, 


0x924000, 


0xA28000, 


0xC30000 ] , 


8 


[ OxDOOSCl, 


0xDlA042, 


0xDlA084, 


0xDlA108, 


0xDlA210, 


0x01A3E0, 


9 


0x00E400, 


0x01C800, 


0x017000, 


0x720000, 


OxE40000, 


0xB80000 ] , 


10 


[ 0x7800El, 


0x88F022, 


0x88F044, 


0x88F088, 


OxFOFOFO, 


0x78E100, 


11 


0x78D200, 


0x78B400, 


0x787800, 


0x990000, 


OxAAOOOO, 


OxCCOOOO ] , 


12 


[ 0xE24031, 


0x738012, 


0x738024, 


0x91C038, 


0x938C40, 


0xElC480, 


13 


0xElC900, 


0x724E00, 


0x020000, 


0x036000, 


0xB40000, 


0xD80000 ] , 


14 


[ 0xCC6009, 


0x66A00A, 


OxAACOOC, 


0xC6C090, 


0x6A60A0, 


OxACAOCO, 


15 


0x6CC900, 


0xA66A00, 


OxCAACOO, 


OxOOFOOO, 


OxOFOOOO, 


OxFOOOOO ] , 


16 


[ OxOOOOBl, 


0x0000E2, 


0x000074, 


0x000008, 


0x7E8100, 


0x828200, 


17 


0x848400, 


0x888800, 


0x909000, 


OxAOAOOO, 


OxCOCOOO, 


OxFFOOOO ] , 


18 


[ OxOOOOBl, 


0x0000E2, 


0x000074, 


0x000008, 


OxOOBlOO, 


0x00E200, 


19 


0x007400, 


0x00D800, 


OxBlOOOO, 


0xE20000, 


0x740000, 


0xD80000 ]] 



20 

21 desd24:= 

22 [LinearCode<GF(2) ,241 [Prune(Intseq(n+0xl000000,2)) : n in code]> 

23 : code in desd24genmats] ; 

The function subcodes takes a doubly even code C containing R as 
an argument, and returns the list of subcodes of codimension 1 of C 
satisfying C D i? up to the action of Aut(C). 

24 subcodes :=funct ion (C,R) 

25 A:=AutomorphismGroup(C) ; 

26 P : =PemutationModule (A , GF (2) ) ; 

27 DC : =Dual (C) ; 

28 DR:=Dual(R); 

29 PDC:=sub<P I VectorSpace(DC)>; 

30 PDR,e:=sub<P | VectorSpace(DR)>; 

31 M,p:=quo<PDR | PDO; 

32 G : =Matr ixGroup (M) ; 

33 X:=[DR| o[l] 00 p e : o in Orbits(G) I not in o] ; 
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34 overcodes:=[sub<DR|DC,x> : x in X] ; 

35 return [Dual(CC) : CC in overcodes] ; 

36 end function; 

Given a sequence of pairs of a code and a number, the function uptoequivalenceDE 
returns a subsequence of complete representatives of codes up to equiv- 
alence with the largest numbers appearing in the second components. 

37 uptoequivalenceDE : =f unction (Ds) 



38 Css : = [] ; 

39 for D in Ds do 

40 ord := #AutomorphismGroup(D[l] ) ; 

41 we := WeightEnuinerator(D[l] ) ; 

42 if not exists (v){i:i in [l..#Css]| 

43 Css [i] [1] eq ord and 

44 Css[i] [2] eq we and 

45 IsEquivalent(Css[i] [3] , D[l])> then 

46 Append(~Css, <ord, we, D[l], D[2]>); 

47 else 

48 Css [v] [4] : =Max ( [Css [v] [4] , D [2] ] ) ; 

49 end if; 

50 end for; 

51 return [<D[3] ,D[4]>: D in Css ]; 



52 end function; 



Basic operation for codes. Given a pair of vectors, the functions 
entrywiseProduct and CstarC return Ci * C2 = Ci fl C2 as the support 
and C * C = (ci * C2 I Ci, C2 e C) respectively. 

53 entrywiseProduct :=func<x,y I 

54 CharacteristicVector(Parent(x) , Support (x) meet Support (y))>; 

55 

56 CstarC :=f unction (D) 

57 k:=Dimension(D) ; 

58 CC:=LinearCode<GF(2) , Length (D) | 

59 [entrjrwiseProduct(D.i,D. j) :i,j in [l..k] | i It j] cat 

60 [D.i : i in [1. .k]]>; 

61 return CC; 

62 end function; 

Given codewords x,y oi a doubly even code C, the functions QForm 
and BForm return Q{x) and B{x,y) respectively. 

63 QForm :=func<u I GF (2) ! (Weight (u) div 4)>; 

64 BForm:=func<u,v|GF(2) ! (#(Support(u) meet Support (v)) div 2)>; 

Given a vector x and a doubly even code D with a basis {ui,U2, • • • , Uk}, 
the function BFormArray returns an array {B{x,Ui))i. Given doubly 
even codes C, D with respective bases {ui, U2, ■ ■ ■ , Uk} and {^1,^2, • • • , vi}, 
the function BFormMatrix returns a matrix (B(ui,Vi))ij. 
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65 BFoirniArray:=function(x,D) 

66 kD:=Dimension(D) ; 

67 return [BFom(x,D. j) : j in [l..kD]]; 

68 end function; 

69 

70 BForTiiMatrix:=function(C,D) 

71 kC:=Dimension(C) ; 

72 kD:=Dimension(D) ; 

73 M:=Matrix(GF(2), kC, kD, 

74 [BFormArrayCC.i, D) : i in [l..kC]]); 

75 return M; 

76 end function; 

Given a doubly even code C, the functions Cmeetrad and CmeetRad 
return the subcode C fl rad C and C fl Rad C respectively, applying 
Lemma [61 

77 Cmeetrad :=funct ion (C) 

78 D:=Dual(CstarC(C)) meet C; 

79 H:=VectorSpace(GF(2) , Dimension (C) ) ; 

80 VD : =VectorSpace (D) ; 

81 g:=hom<VD->H|BFormMatrix(D,C)>; 

82 rad : =sub<D | Kernel (g) > ; 

83 return rad; 

84 end function; 
85 

86 CmeetRad :=funct ion (C) 

87 rad:=Cmeetrad(C) ; 

88 k:=Dimension(rad) ; 

89 H:=VectorSpace(GF(2) ,1) ; 

90 VD : =VectorSpace (rad) ; 

91 g:=hom<VD->H| [ [QForm(rad. i)] : i in [l..k]]>; 

92 Rad : =sub<rad | Kernel (g) > ; 

93 return Rad; 

94 end function; 

Doubly even codes which contain each triply even radical. The 

function outsideVectors returns a complete list of representatives of 
cosets (C * * C)-^ n C) up to the action of Aut(C). Given a 

doubly even code C, the function existsOutsideRad returns true if 
and only if Rad C C, applying Lemma [23] and Lemma [251 

95 outsideVectors :=function(C) 



96 


U:=Generic(C) ; 


97 


A:=AutomorphismGroup(C) ; 


98 


P : =PemutationModule (A , GF (2) ) 


99 


D : =Dual (CstarC (C) ) ; 


100 


E:=C meet D; 


101 


PD,e:=sub<P | VectorSpace(D)> 


102 


PDC:=sub<P 1 VectorSpace(E)>; 


103 


M,p:=quo<PD | PDO; 


104 


G : =Matr ixGroup (M) ; 
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105 return {D!(o[l] 00 p @ e) : o in Orbits (G) I not in o>; 

106 end function; 

107 

108 existsOutsideRad:=function(C) 

109 H:=VectorSpace(GF(2) , Dimension (C) ) ; 

110 D:=Dual(CstarC(C)) meet C; 

111 VD : =VectorSpace (D) ; 

112 g:=hom<VD->H| BFormMatrix(VD,C)>; 

113 Im:=Image(g) ; 

114 rad := Kernel (g); 

115 bl :=exists(u){ i : i in [1 . .Dimension(rad)] I QFormCrad. i) ne 0}; 

116 X : =outsideVectors (C) ; 

117 b2:=exists(v){x : x in X | 

118 imgx in Im and (bl or QForm(x+imgx 00 g) eq 0) 

119 where imgx:= H!BFormArray(x,C)}; 

120 return b2; 

121 end function; 

The record RF equips the following objects for a doubly even code of 
length 24. 



122 


RF : =recf omat< 


123 


c, 


// the original code 


124 


R, 


// the triply even radical of C 


125 


prd. 


// the max dim of radical of supcode of codim = 1 


126 


CR, 


// the quotient space C/R 


127 


p. 


// the projection C -> C/R 


128 


X, 


// the array [ x in CR 1 Q(x) = ] 


129 


px. 


// the projection V(X)->C/R 


130 


CC, 


// the triply even check code 


131 


AutCR 


// Aut(C) meet Aut(R) 


132 


>; 





Given a doubly even code C and its triply even radical i?, the procedure 
profiles constructs the quotient C/R, the projection p : C ^ C/R, 
the singular points X, the automorphism group Aut(C) fl Aut(i?) and 
the triply even check code, and then returns a record containing them. 

133 prof iles:=f unction (C, prd) 

134 s:=rec<RF | C:=C, prd:=prd>; 

135 s ' R : =CmeetRad (C) ; 

136 s'CR,s'p:=VectorSpace(C)/VectorSpace(s'R) ; 

137 s'X:=[x:x in s'CR|QForm(x 00 s'p) eq 0]; 

138 M:=Matrix(GF(2), #s'X, Dimension(s'CR) , s'X); 

139 s ' px : =hom<VectorSpace (GF (2) , #s ' X) ->s ' CR I M> ; 

140 s' AutCR :=AutomorphismGroup(C) 

141 meet AutomorphismGroup(s'R) ; 

142 s ' CC : =LinearCode (Kernel ( s ' px) ) ; 

143 return s; 

144 end function; 

The procedure constAllSubcodeContainsRad constructs the list of all 
doubly even codes of length 24 containing its triply even radical. 
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145 constAllSubcodeContainsRad: =f unction (niaxcodes24) 



146 codes :=[[ profiles (D, 0) : D in maxcodes24] ] ; 

147 print "=> Now, constmcting all admissible doubly even codes of length 

24. . ."; 

148 for i in [1. .9] do 

149 d:=12-i; 

150 ovcodes : =codes [#codes] ; 

151 reps:=&cat[[<C,Dimension(s'R)>:C in subcodes(s'C, s'R)] : s in ovcodes]; 

152 reps :=uptoequivalenceDE (reps) ; 

153 reps:=[S : S in reps | not existsDutsideRad(S[l] )] ; 

154 printf "=> Completed for dim=yo3o, the number of codes=%4o . \n" , d, #reps 

155 Append(~codes, [prof iles(D[l] , D[2]) :D in reps]); 

156 if IsEmpty(reps) then 

157 break i; 

158 end if; 

159 end for; 

160 printf "=> This is the expected result : XoAn", 

161 [#x:x in codes] eq [9,42,160,377,437,220,36,1, 0]; 

162 return &cat (codes) ; 



163 end function; 

Identification of mctximal triply even codes. Given a triply even 
code, the function isMcLximal returns true if and only if the code is a 
maximal triply even code. 

164 isMaximal:=function(C) 



165 D:=CodeComplement(Dual(CstarC(C)), C) ; 

166 t : =exists (u) {x : X in D| x ne and QForm(x) eq 

167 and forall(v){i:i in [1 . .Dimension(C)] I BForm(x,C.i) eq 0»; 

168 return not t; 



169 end function; 

The procedure appendCode appends a new maximal triply even code 
to the list of codes. 

170 appendCode :=procedure(~codenum, "maxcodes, reps, Ds, id) 



171 codenum:=codenum + #Ds; 

172 for D in Ds do 

173 if isMaximal(D) then 

174 Append(~maxcodes, D) ; 

175 invt : =<Dimension (D) , NumberOf Words (D , 8) > ; 

176 idO : =Position(reps [2] , invt) ; 

177 if idO ne and not IsEquivalent(D,reps[l] [idO]) then 

178 idO := 0; 

179 end if; 

180 printf 

181 "Foimd a MTE code = Rep.%2o of dim=°/.2o from DE code No.°/.o : %o.\n", 

182 idO, Dimension(D) , id, reps [3] [idO+1] ; 

183 end if; 

184 end for; 



185 end procedure; 
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Triply even codes constructed from the combinations with 
Ci = C*2- We enumerate all codes obtained from the method in Propo- 
sition [19] with Ci = C2. 

Given a doubly even code and its triply even radical, the function 
constDoubleCosetsCC returns the representatives of double cosets 

(33) g,{c,R)\g,{c,R)/go{c,R). 

186 constDoubleCosetsCC :=function(s) 

187 CRs:={@x:x in s'CRix ne 00}; 

188 SCRs:=Sym(CRs) ; 

189 GO:=sub<SCRs|{[((x M s'p)"g) s'p:x in CRs] : 

190 g in Generators (s'AutCR)}>; 

191 GLCR:=sub<SCRs|{[x~g:x in CRs]: 

192 g in Generators (GL(s'CR))>>; 

193 Gl:=Stabilizer(GLCR, {x : x in s'X | x ne 0}); 

194 return DoubleCosetRepresentatives(Gl, GO, GO); 

195 end function; 

Given a doubly even code C and the double cosets fl33l) . the function 
result ingCodesCC returns triply even codes constructed from the code 
C using the method in Proposition [12] with Ci = €2 = C. 

196 resultingCodesCC:=function(s, dc) 

197 D:=DirectSum(s'R, s'R); 

198 k:=Dimension(s'CR) ; 

199 Ml:=Matrix( [s'CR.i @@ s'p : i in [l..k]]); 

200 codes :=[D+LinearCode(HorizontalJoin (Ml, M2)) where 

201 M2:=Matrix([((s'CR.i)"g) @@ s'p : i in [l..k]]) 

202 : g in dc] ; 

203 return codes; 

204 end function; 

The object partsDB is the set of doubly even codes of length 24 con- 
taining its triply even radical. The function duplextype returns the 
list of all maximal triply even codes and the number of triply even 
codes of length 48 constructed from partsDB with 6*1 = 6*2. 

205 duplextype :=f unction (partsDB, reps) 



206 maxcodes : = [] ; 

207 codenum : =0 ; 

208 excodenum : =0 ; 

209 for id in [1 . .#partsDB] do 

210 s:=partsDB[id] ; 

211 k:=Dimension(s'CR) ; 

212 if k eq then 

213 if Diniension(s'R) eq s'prd then 

214 excodenum :=excodenum+l; 

215 else 

216 D:=DirectSum(s'R, s'R); 

217 appendCode( "codenum, ~maxcodes, reps, [D] , id) ; 

218 end if ; 
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219 else 

220 doubleCosets : =constDoubleCosetsCC(s) ; 

221 if Diniension(s'R) eq s'prd then 

222 Remove ("doubleCosets, 1); 

223 excodenum:=excodenum+l; 

224 end if ; 

225 if not IsEmpty (doubleCosets) then 

226 list:=resultingCodesCC(s, doubleCosets); 

227 appendCode(~codenum, ~maxcodes, reps, list, id); 

228 end if; 

229 end if; 

230 end for; 

231 return maxcodes, codenum, excodenum; 



232 end function; 

Triply even codes constructed from the combinations with 
Ci ^ (72. We enumerate all codes obtained from the method in Propo- 
sition [Tl with Ci ^ C2. 

Given a pair of doubly even codes and an isomorphism between their 
triply even check codes, the function isometry returns an isometry 
between them. 

233 isometry:=function(sl, s2, g) 



234 CX:=Image(sl'px) ; 

235 BCX, p:=quo<sl'CR|CX>; 

236 bCRl:=Basis(CX) cat [x 00 p : x in Basis (BCX) ] ; 

237 bCX2 : = [bCRl [i] 00 sl'px g s2'px : i in [1. .Dimension(CX)]] ; 

238 BbCX2, p:=quo<s2'CR|bCX2>; 

239 bCR2:=bCX2 cat [x 00 p : x in Basis (BbCX2) ] ; 

240 return hom<sl'CR->s2'CR | [bCRl [i] ->bCR2 [i] : i in [1. .#bCRl]]>; 



241 end function; 

Given the object partsDB, which is the set of doubly even codes of 
length 24 containing its triply even radical, the function isometricPairsClC2 
returns the list of isometric pairs of distinct doubly even codes and an 
isometry between them. 

242 isoinetricPairsClC2 : =f unction(ss) 



243 ClC2s:=&cat[[<i, j, isometry(ss[i] ,ss[j] ,g)> 

244 : j in [i+1 • • #ss] 

245 I Dimension(ss [i] 'CR) eq Diniension(ss [j] 'CR) 

246 and #ss[i]'X eq #ss[j]'X and isEq 

247 where isEq, g := IsEquivalent(ss [i] 'CC, ss[j]'CC)] 

248 : i in [1. .#ss]] ; 

249 printf "The number of hybrid pairs = 125: °/oO.\n", #ClC2s eq 125; 

250 return ClC2s; 



251 end function; 

Given a pair of doubly even codes and an isometry, the function constDoubleCosets 
returns the double cosets 



(34) h-'goiC2, R2)h\G,{C^, Ri)/Qo{Ci, Ri). 
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252 constDoubleCosetsClC2:=function(sl, s2, h) 



253 CRs:={@ X : X in sl'CR | x ne @>; 

254 SCRs:=Sym(CRs) ; 

255 G01:=sub<SCRs | {[((x 00 sl'p)"g) sl'p 

256 : X in CRs] : g in Generators (sl'AutCR)}>; 

257 002 : =sub<SCRs | {[((x h 00 s2'p)"g) s2'p 00 h 

258 : X in CRs] : g in Generators (s2'AutCR)}>; 

259 GLCR:=sub<SCRs | {[x"g : x in CRs] : g in Generators (GL(sl'CR))>>; 

260 Gl:=Stabilizer(GLCR, {x : x in sl'X | x ne 0}); 

261 return DoubleCosetRepresentatives(Gl, GOl, G02) ; 



262 end function; 

Given a pair of doubly even codes Ci and C2, an isometry h from Ci/Ri 
to C2/R2 and the double cosets flMl) . the function resultingCodesClC2 
returns triply even codes constructed from the pair of codes using the 
method in Proposition [191 

263 resultingCodesClC2:=function(sl, s2, h, dc) 



264 k:=Dimension(sl'CR) ; 

265 D:=DirectSum(sl'R, s2'R) ; 

266 Ml : =Matr ix ( [sl'CR. i 00 sl'p : i in [l..k]]); 

267 codes :=[D+LinearCode(HorizontalJoin (Ml, M2)) where 

268 M2:=Matrix([((sl'CR.i)''g) h 00 s2'p : i in [l..k]]) 

269 : g in dc] ; 

270 return codes; 



271 end function; 

Recall that the object partsDB is the set of doubly even codes of length 
24 containing its triply even radical. The function hybridtype returns 
the list of all maximal triply even codes and number of triply even 
codes of length 48 constructed from partsDB with Ci ^ C^- 

272 hybridtype :=f unction (partsDB, reps) 



273 maxcodes : = [] ; 

274 codenum : =0 ; 

275 clc2s:=isometricPairsClC2 (partsDB) ; 

276 for id in clc2s do 

277 sl : =partsDB [id [1] ] ; 

278 s2 : =partsDB [id [2] ] ; 

279 h:=id[3]; 

280 k:=Dimension(sl'CR) ; 

281 if k eq then 

282 D:=DirectSum(sl'R, s2'R); 

283 appendCode(~ codenum, "maxcodes, reps, [D] , <id[l] ,id[2]>) ; 

284 else 

285 doubleCosets:=constDoubleCosetsClC2(sl, s2, h) ; 

286 list:=resultingCodesClC2(sl, s2, h, doubleCosets) ; 

287 appendCode(~ codenum, "maxcodes, reps, list, <id[l] , id[2]>); 

288 end if; 

289 end for; 

290 return maxcodes, codenum; 



291 end function; 
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Representative examples of mciximal triply even codes. We 
give 10 maximal triply even codes {I'(C) | C € A} and C'(Tio) of 
length 48. These codes are constructed by the functions tildeD and 
TriangularGraphCode. 

292 tildeD := fimction(C) 

293 R := CmeetRad(C) ; 

294 return Juxtaposition(C,C)+DirectSiim(R,R) ; 

295 end function; 

296 

297 TriangularGraph:=function(v) 

298 X:=SetToIndexedSet (Subsets ({1 . .v>,2)) ; 

299 return #X, Matrix (GF (2) , #X, #X, [[#(x meet y) :y in X] : x in X]); 

300 end function; 

301 

302 TriangularGraphCode :=f unction (v) 

303 n, M:=TriangularGraph(v) ; 

304 r:=(-n) mod 8; 

305 return PadCode(LinearCode(M) ,r) + RepetitionCode(GF(2) , n+r) ; 

306 end fimction; 

The object repMTECodes is the hst of 10 maximal triply even codes 
equipped with their dimensions and the numbers of their codewords of 
weight 8. 

307 repMTECodesl:=[ tildeD(C) : C in desd24 ] cat [TriangularGraphCode (10) ] ; 

308 repMTECodes2:=[<Dimension(C) ,NumberDfWords(C,8)> : C in repMTECodesl] ; 

309 repMTECodes3: = ["New!", "tD( g_{24} )", "tD( d_{24}'-{+} )", 



310 "tD( d_{12>-{2+> )", "tD( (d_{10>e_7"2)-{+} )", 

311 "tD( d_8"{3+} )", "tD( d_6-{4+} )", "tD( d_4-{6+} )", 

312 "tD( d_{16>~{+}\oplus e_8 )", "tD( e_8"{\oplus3}\} )", 

313 "tT_{10>"] ; 

314 dim_repMTECodes:={* Dimension (C) : C in repMTECodesl*}; 

315 printf "Representative codes are inequivalent each other: %o.\n", 

316 #repMTECodesl eq #Seqset (repMTECodes2) and 

317 dim_repMTECodes eq {* g^'l, 13~"7, 14""1, 15""! *}; 



318 repMTECodes : =<repMTECodesl , repMTECodes2 , repMTECodes3> ; 

Non existence of the other mcLximal triply even code. In this 

subsection, we aim to ensure that there does not exist any maximal 
triply even code of length 48 except for the representative examples in 
the previous subsection up to equivalence. 

First, we enumerate all doubly even codes of length 24 which contain 
their triply even radicals. 

319 partsDB:=constAllSubcodeContainsRad(desd24) ; 

320 table:= [ [Integers !0: j in [1..13-k]]:k in [1..9]]; 

321 for k in [1 . .#partsDB] do 

322 i : =13-Dimension(partsDB [k] ' C) ; 

323 j : =Dimension(partsDB [k] ' R) ; 

324 table[i] [j] + :=l; 
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325 end for; 

326 printf "The number of acJmissible codes is same as expected: °/oO.\n", 

327 table eq 

328 [ 

329 [ 7, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, ], 

330 [ 33, 6, 3, 0, 0, 0, 0, 0, 0, 0, ], 

331 [ 130, 19, 10, 1, 0, 0, 0, 0, 0, ], 

332 [ 308, 40, 23, 5, 0, 1, 0, 0, ], 

333 [ 363, 37, 25, 10, 1, 1, 0, ], 

334 [ 180, 16, 10, 11, 2, 1, ] , 

335 [ 27, 2, 0, 4, 2, 1 ] , 

336 [ 0, 0, 0, 0, 1 ] , 

337 [ 0, 0, 0, ] 

338 ] ; 

Second, we check the maximahty of triply even codes constructed from 
all the doubly even codes in duplicate. 

339 duplex_max, duplex_num, exduplex_num 

340 :=duplextype(partsDB, repMTECodes) ; 

341 printf "°/o3o maximal codes of duplex type found. \n",#duplex_max; 

342 printf "This is the expected result : 7oO . \n" , 

343 <#duplex_max, duplex_num, exduplex_num> eq <30,214,1268>; 

Next, wc check the maximality of triply even codes constructed from 
all the pairs of distinct doubly even codes. 

344 hybrid_max, hybrid_num:=hybridtype(partsDB, repMTECodes); 

345 printf "7o3o maximal codes of hybrid type found. \n" ,#hybrid_max; 

346 printf "This is the expected result: "/oO.Xn", 

347 <#hybrid_max, hybrid_nuiii> eq <5,225>; 

Result. A classification of triply even codes of length 48 has been 
completed into the 10 codes. This calculation has been completed in 
the total time: 650.240 seconds, the total memory usage: 534.91MB 
under the environment using "Intel® Core^M 2 Duo CPU T7500 @ 
2.20GHz". 
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